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ABSTRACT 

We present the results of both analytical and numerical calculations of the amplitude 
of the reflection component in X-ray spectra of galactic black hole systems. We take 
into account the anisotropy of Compton scattering and the systematic relativistic bulk 
motion of the hot plasma. In the case of single scattering approximation the reflection 
from the disc surface is significantly enhanced due to the anisotropy of Compton 
scattering. On the other hand the calculations of multiple scattering obtained using 
the Monte Carlo method show that the anisotropy effect is much weaker in that case. 
Therefore, the enhanced backscattered flux may affect the observed spectra only if 
the disc surface is highly ionized, which reduces the absorption in the energy band 
corresponding to the first Compton scattering. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - X-rays: 
galaxies - X-rays: stars 
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1 INTRODUCTION 

Hard X-ray spectra of the galactic black hole systems are 
well described by a power law primary emission along with 
the pronounced reflected component, which causes the ob- 
served flattening of the spectrum. The primary emission is 
likely to be produced by Compton upscattering of soft pho- 
tons on thermal electrons in hot, optically thin medium close 
to a relatively cold accretion disc, being the source of seed 
photons for Comptonization (see e.g. review in Poutanen 
1998). A fraction of the upscattered photons is directed to- 
wards the disk and can be reflected from its surface, giv- 
ing the rise to reflected continuum and fluorescent iron line 
emission (Lightman & White 1988; George & Fabian 1991). 

Observational data for Cyg X-1 and other black hole 
systems in their hard/low state show often rather hard spec- 
tra (photon spectral index F ~ 1.5 — 1.9; Poutanen et al. 
1997; Gierlihski et al. 1997; Dove et al. 1997; see also Pouta- 
nen 1998), while the amplitude of reflection R covers the 
broad range of values between and 2. Moreover, R and F 
are correlated (Zdziarski, Lubihski & Smith 1999; Revnivt- 
sev, Gilfanov & Churazov 1999), in the sense that the harder 
the spectrum, the smaller the amplitude of reflection. The 
correlation exists both within the low/hard state and when 
sources change their spectral state (Zycki, Done & Smith 
1998). 

These observations cannot be explained by the model 
in which static, continuous corona covers the cold accretion 
disc, as it predicts the power law slope F ^ 2 (i.e. rather soft 



spectra) and the reflection amplitude _R ~ 1.0. Among the 
possible models, which could reproduce the reflection am- 
plitude in the range R = — 1 there are two competitive: (i) 
cold disc disrupted in the inner part (e.g. Poutanen, Krolik 
& Ryde 1997; Esin, McClintock & Narayan 1997) and (ii) 
a highly ionized, non-disrupted disc (Nayakshin, Kazanas & 
Kallman 2000; Ross, Fabian & Young 1999). Detailed shape 
of the reflected continuum depends on the geometry, ion- 
ization state and abundances of elements in the scattering 
medium. However, spectral fitting does not always allow to 
constrain these parameters independently. It is possible to 
explain the observed spectra of GBH in terms of weakly ion- 
ized, or neutral reflection from the disc which inner radius 
is of the order of 50 Rg (e.g. Done & Zycki 1999 for the hard 
state of Cyg X-1), as well as with highly ionized reflection 
from the disc extending to the marginally stable orbit (as 
suggested by Ross et al. 1999 for the hard state of Cyg X-1; 
see also Done & Nayakshin 2000). 

The third possible model, in which both the reflection 
amplitudes R > 1 and R < 1 are possible, is a mildly rel- 
ativistic outflow/inflow in the corona (Beloborodov 1999). 
Relativistic aberration reduces the hard X-ray flux scattered 
towards the disc, which leads to reduction of the reflected 
component and the soft flux from reprocessing entering the 
corona. In order to obtain quantitative agreement with ob- 
served spectral indices and reflection amplitudes, the model 
requires additional reduction of the soft flux intercepted by 
the hot plasma. This leads to the 'active regions' geome- 
try (e.g. magnetic flares above the disc; Haardt, Maraschi & 
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Ghisellini 1994) rather than a continuous corona. However, 
the coronal outflow may reproduce the reflection amplitude 
R < 1 and in order to obtain i? > 1 an inflow of the plasma 
must be postulated. 

In this article we reanalyze the non-static corona model, 
with plasma moving at relativistic speed in the direction per- 
pendicular to the disc surface. We study the dependence of 
the amplitude of the reflection component on the bulk mo- 
tion velocity, taking into account the thermal motion of elec- 
trons within the plasma determined by the electron temper- 
ature Te and we discuss the resulting shape of the spectra. 
We also check if high reflection amplitudes {R > 1) might 
be explained in the frame of the outflow model but tak- 
ing into account possible high ionization of the disc surface 
and the anisotropy of Compton scattering. We emphasize 
the importance of the first Compton scattering, which plays 
crucial role at lower energies. For highly ionized disc surface 
the effect of absorption by heavy elements is reduced and 
the reflected spectrum makes a substantial contribution to 
the total spectrum in the ~ 0.5 — 5 keV band. Therefore 
we present firstly the semi-analytical calculations in the ap- 
proximation of single scattering, and after that we perform 
numerical simulations of multiple scattering, which is re- 
sponsible for the power law shape of the hard X-ray tail. 

The contents of the paper is as following. In Section ti 
we analyze semi-analytically the amplitude of the refiection 
in a single scattering approximation (after Ghisellini et al. 
1991) but taking into account anisotropy of scattering within 
the rest frame of the electron, describing the thermal motion 
without the assumption of highly relativistic beaming and 
incorporating the systematic bulk motion (outfiow) of the 
corona. Since the effect of multiple scattering is essential in a 
real situation, in Section H we repeat the computations for a 
corona of a given optical depth and an electron temperature 
using a Monte Carlo method for both outwards and inwards 
directions of bulk velocity. The dependence of the refiec tion 
amplitude on the outflow velocity is presented in Section 3.1 
The spectra resulting fro m Monte Carlo computations are 
presented in Section 3.2, The discussion of the results is 
given in Section m 



2 SINGLE SCATTERING APPROXIMATION 
FOR THE AMPLITUDE OF REFLECTION 

In this section we generalize the determination of the ampli- 
tude of the reflection component derived by Ghisellini et al. 
(1991) for a slab geometry of a hot corona and anisotropic 
soft photon input from underlying disc. Those results were 
obtained assuming isotropic electron scattering in the rest 
frame of an electron and relativistic chaotic motion of elec- 
trons. We introduce subseque ntly the anisotropy of the 
Thomson cross-section (Section 2.1), we relax the assump- 
tion tha t th e thermal motion of the plasma is relativistic 
(Section 2.2), and flnally we introduce the systematic bulk 



motion of the corona (Section 

We show, that including the effect of anisotropic soft 
photon distribution in Compton scattering process the re- 
flection amplitude R ^ 2 can be obtained. It's value is 
somewhat larger when the angular dependence of Thom- 
son cross-section is taken into account. On the other hand 
for low electron velocities (7 < 2) the angular distribution of 



scattered radiation becomes important and the value of R is 
reduced. Adding the systematic coronal outflow reduces the 
anisotropy effect only in the range of high bulk velocities. 
However all these results stand mostly for the flrst scatter- 
ing, while for the multiple scattering in the cloud of optical 
depth r ~ 1 the reflection is weaker (see Section H). 



2.1 The effect of angular dependence of the 
Thomson cross-section 

In this section we calculate the power of Compton radia- 
tion scattered by electrons with isotropic relativistic veloc- 
ity held V = f3c, /3 ~ 1. This approach is appropriate for a 
non-thermal plasma. The radiation held is anisotropic and 
incoming photons subtend a restricted solid angle. We as- 
sume the geometry is the same as in Ghisellini et al. (1991). 
However, in our calculations we take into account the angu- 
lar dependence of Thomson cross-section (Rybicki & Light- 
man, 1979). Therefore the emitted power is given by 



n«,7) = ^ 



Pi. 



8(1 + /3V3) 



27r — 0„i„ 



{1- f3 COS 9 f{l + cos^ 9) sin 9d9 



(1) 



Here Piso is the power of isotropic emission and P{a) is 
the power of radiation scattered by electron of the angle a 
between its velocity and the axis of symmetry of incoming 
photons. The limits ^min, ^max and (^minare functions of a 
and are given by equations (see Ghisellini et al., 1991): 



max(0, a — vr/2) 
: min(a 4- 7r/2, n) 



arccos[(tan 9 tan a) 



(2) 



for < e < 7r/2 - Q 
or 37r/2 - a < < tt 
otherwise 



Here we follow the assumption of Ghisellini et. al (1991) , 
that all photons are emitted in the direction of electron's 
velocity. This is justified in the case of 7 3> 1, when the 
direction of the motion of an electron (i.e the angle a) can 
be identified with the viewing angle. In Figure hi we plot 
the ratio of the total power emitted by electrons moving 
downwards to that of electrons moving upwards: 



P(7) 



/J/2^(")^i"(")'^" 



"J. 

ir/2 



[^ P{a) sin(Q:)dQ 



(3) 



It is worth seeing that this ratio saturates for large 7 at 
the value of i? ~ 2.55, which is somewhat larger than that 
obtained under assumption of da/dO. = gt/^-k {R ~ 2.2). 
This result might be of possible importance for non-thermal 
Comptonization models. It seems however to be in conflict 
with the recent results of spectral fitting of the data for the 
soft state of Cyg X-1 (Gierliiiski et al. 1999). 

In Figure y we plot the ratio of the power emitted down- 
wards to that emitted by electron moving in one particular 
direction, given by 
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Figure 1. The ratio of tlie total power emitted by electrons mov- 
ing downwards to that of moving upwards as a function of elec- 
tron's Lorentz factor 7. 
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Figure 2. The ratio of the total power emitted by electrons mov- 
ing downwards to that of electron moving in one particular direc- 
tion, as a function of electron temperature kT[keV], for an incli- 
nation angle a = 30° , under the assumption of efficient beaming 
(dashed line), and for the same inclination angle (i = 30°) but av- 
eraged over the whole range of electron's velocity direction (solid 
hne). 



R{a,kT) = 



/J/2 ^^^^ sm{a)da 



P{a) 



(4) 



This can be identified with the reflection amphtude seen 
by the observer with the viewing angle a. We plot the re- 
sult versus the electron temperature instead of 7 in order 
to expose better the lower velocity part of the curve but 
we remind that those results (dashed line) are not supposed 
to be accurate for not very highly relativistic velocities. We 
see that for low inclination angle the expected enhancement 
of the reflection is very large, up to a factor 3.5 for typical 
plasma temperatures. 



2.2 The effect of mildly relativistic thermal 
motion 

If the plasma is rather thermal than non-thermal, the elec- 
tron velocity is only mildly relativistic and the assumption 
of all the emission being beamed in the direction of motion 
of an electron used in Section 2.1 is not valid. Therefore, 



in this section we consider the angular distribution of the 
scattered radiation, given by 



P(Q,7,9out) = 



P(a,7) 



274(1 -h/3coseout)3 



(5) 



where Oout ~ tt means scattering in the direction of elec- 
tron's movement. The accretion disc receives a fraction of 
scattered radiation that depends on the direction and veloc- 
ity of electron 



PdiBc{a,j) 



©" 



27r — 0" 



P(a, 7, 6out) sin 9outd9out#out (6) 



where OoS'', ©ouT and (?iou" are given by equations (g). 

The viewing angle i in this case cannot be identified 
with the direction of electron motion given by a. Instead, 
we have a relation 



cos 9ou 



(sin i sin ip sin a + cos i cos a) 



(7) 



Here a and ip determine the electron's velocity vector direc- 
tion: 



ii(cos ip sin a, sin ip sin a, cos a) 



(8) 



In Figure H we show the assumed geometry scheme. Note, 
that the angle 0out is measured from the vector —v to the 
direction of emitted photon, and therefore we keep the in- 
tegration limits given by equations (B) unchanged (see also 
Fig. 1 and Fig. 2 in Ghisellini et al. 1991). 

We calculate the reflection amplitude averaged over the 
whole range of a and ip as & function of electron velocity 
and viewing angle, according to the formula 



7?(7,i) = 



1 Jq Jq " Pdiac{a, 7) sin adadip 

2^ Jo" /o " Pi'^' 7. ©out (a, ifi, i)) sin adadip 



(9) 



In Figure W we plot the refiection amplitude as a function of 
7 for three different values of viewing angle. The plot for an 
inclination i = 60° roughly corresponds to the value aver- 
aged over all inclinations, as presented in Figure hi for non- 
thermal plasma. We see that the anisotropy is only slightly 
reduced if the total coUimation assumption is relaxed. This 
effect is mostly seen in Figure where we present with the 
solid curve low velocity (moderate temperature) tail of the 
distribution. 



2.3 Mildly relativistic bulk motion 

In this section we assume that bulk velocity vector is perpen- 
dicular to the disc surface and directed outwards. We calcu- 
late the net electron velocity as a sum of thermal chaotic mo- 
tion and systematic (bulk) outflow. The angle between the 
net velocity vector and vertical axis, a , is connected with 
angle a via relativistic velocity transformation. Therefore 
the net Lorentz factor, 7', depends on the angle a' as well 
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Figure 3. The assumed geometry scheme. The electron is in the 
center of the reference frame and has the velocity v at the in- 
clination a to the symmetry axis of the disc and at the angle ip 
measured in the disc plane. The photon is emitted at the angle 
i to the disc axis, and at the angle (©out — i") to the electron's 
velocity vector. 
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Figure 4. The reflection amplitude for the viewing angle i = 
(solid line), 30° (short-dashed line) and 60° (long-dashed line) as 
a function of electron's Lorentz factor 7 in case of pure thermal 
motion. 



as on bulk and thermal velocities. In this case the amount 
of reflection is given by: 
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Figure 5. The reflection amplitude for the viewing angle i = 
0° (boxes) and i = 45° (triangles) as a function of /3buik for 
electron temperature Ahorm = 0.63 (top panel) and /Sthcrm = 
0.72 (bottom panel). 



viewing angle i. The electron temperature corresponds to a 
single value of electron velocity, as for velocity distribution 
the analytic approximation would not work. 

In the case of Vtherm = we obtain the same solution 
as in Beloborodov (1999). In Figure owe plot this solution, 
calculated for different viewing angles. 

The comparison of Figures n and n shows that single 
scattering approach predicts very strong dependence of the 
amplitude of reflection on the plasma temperature due to 
anisotropy of the Compton scattering. When the thermal 
motions are important (iithcrm/c ~ /3buik) the first reflection 
is significantly enhanced, by a factor of a two, and the higher 
the temperature the stronger the effect. The values of the 
amplitude of reflection cover the whole range between and 
~ 2 for outflow solutions (/3buik > 0). 

The effect of reflection enhancement drops rapidly when 
the bulk velocity is dominant (/3buik > ^'thcrm/c). In this 
case the maximum angle a' is much smaller than n and the 
integration limits in the equation ( y_0|) must be changed. 
The reflection amplitude increases then with the value of 
assumed viewing angle i, while in the case of more significant 
thermal motion the trend is opposite. 



i?(/3bulk, Ahcrm, J 
r2w 



1 

2^ 



/o /o -Pdisc(Q',/3buik, Ahcrm)sina'da'd<^' 

lo /o "" -P("'' /^buik, Ahorm, 0out («', «)) siu a' da' dip' 

In Figure |5| we plot the dependence of the amount of 
reflection on the bulk velocity to the light velocity ratio, 
/3buik, for different values of electron temperature /Sthcrm and 



3 REFLECTED SPECTRA FROM MONTE 
CARLO SIMULATIONS 

The roughly power law shape of the primary emission com- 
ponent in X-ray spectra in GBH is most likely due to the 
effect of multiple scatterings within the hot plasma. 

It is well known from analytical solutions and Monte 
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Figure 6. The reflection amplitude for the viewing angle i = 0° 
(solid line), 30° (short-dashed line) and 60° (long-dashed line) as 
a function of /3bulk for pure bulk motion. 



Carlo simulations that the contribution of the first scat- 
tering to the total spectrum is rather specific (see Stern 
et al. 1995, Svensson 1996, Haardt, Maraschi & Ghisellini 
1997). It means, that when the hard X-ray spectrum forms 
via multiple scatterings, only the first one is infiuenced by 
the anisotropy of seed photon distribution. This is the rea- 
son why any anisotropy effects are present in the low en- 
ergy part of the spectrum. In the thermal medium, even for 
small optical depths, the power law spectrum is shaped by 
multiple scatterings. Therefore the semi-analytical compu- 
tations, dealing with the first scattering process, in the case 
of both thermal and systematic bulk motion of the electrons 
within the corona can only serve as a guide and a help to 
understand the numerical results. The fijlly reliable answer 
can only be provided by full Monte Carlo simulations of the 
Comptonization process within the corona. 

In this Section we compute the amplitude of the Comp- 
ton reflected component using a Monte Carlo comptoniza- 
tion code and we show the corresponding spectra. We con- 
centrate on highly ionized reflector, as the approximations 
used in Beloborodov (1999) may not be valid in that case. 
The code employs standard algorithms for simulating the 
inverse Compton scattering, and it was written following 
descriptions by Pozdnyakov, Sobol & Sunyaev (1983) and 
Gorecki & Wilczewski (1984). Modifications to the code nec- 
essary to implement the bulk motion are described in Ap- 
pendix A. Following Beloborodov (1999) we assume that the 
comptonizing region as a whole is stationary and its geom- 
etry is that of a slab. 



3.1 Reflection amplitude 

The bulk velocity vector is assumed perpendicular to the 
plane of the disc but the direction of the velocity can be 
both outwards and towards the disc. The disc is a source 
of soft photons for the comptonization. Photons backscat- 
tered from the cloud form radiation illuminating the disc. 
The usual Compton reflection process is then simulated by 
another Monte Carlo routine (Zycki & Czerny 1994), as- 
suming the abundances given in Morrison & McCammon 
(1983). The opacities were computed using the code de- 



2.5 



1.5 



0.5 



\ 



-0.5 



0.5 



Figure 7. Amplitude of the Compton-reflected component as 
a function of the bulk outflow velocity, /3bulki computed using 
the Monte Carlo method for multiple scattering in the plasma 
of optical depth Tos = 0.8 and electron temperature kTe = 100 
keV. Open triangles indicate the results for the inclination angle 
cos 4 = 0.3, stars mark the results for cos i = 0.9, and solid squares 
result from the angle averaged spectrum. 



scribed in Done et al. (1992) for the ionization parameter 
^ = Fx/in^r'^) — 10''. Further scattering of the reflected 
photons in the hot comptonizing cloud is not considered, 
since the intercepted fraction would be geometry-dependent 
(e.g. factor ^^ in Beloborodov 1999). 

The refiection amplitu de R is defined here (cf. Be- 
loborodov 1999 and Section 2.3), as the ratio of the energy 
integrated fluxes: 

_ fback (/3bulk) 



R{f3b 



(11) 



.Fdir(/3bulk,«) 

Here F^ir is the flux directly escaping from the comptonizing 
cloud towards an observer, at the inclination angle i, and 
^back is the flux directed towards the disc. We note that 
our definition of R would only be equivalent to that used in 
spectral fitting, if the shapes of spectra used in models fit 
to the data were exactly the same as in our simulations (for 
both the primary and the reflected components). Since this is 
not necessarily the case in practice, translating an amplitude 
R inferred from spectral modelling to a bulk velocity may 
not be accurate. 

The parameters of the comptonizing cloud: the electron 
temperature kTc and optical depth Tos are chosen so that the 
resulting comptonized spectra correspond to typical spectra 
of GBH. Figure M shows the resulting amplitude of the repro- 
cessed component as a function of the plasma bulk velocity 

/3bulk- 

Comparison of the Figures m and a shows that the en- 
hancement of the reflection due to the anisotropy of the 
first scattering is sharply reduced if the subsequent scatter- 
ings are taken into account. Only for rather high inclination 
angles is the obtained reflection amplitude larger than 1.0 
for /3buik > 0, as the larger the viewing angle, the fewer pho- 
tons are able to escape from the slab towards the observer. 
This means, in particular, that bulk velocities directed to- 
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Figure 8. Coinptonizcd spectra resulting from the outflowing 
corona model obtained by Monte Carlo simulations. Parameters 
were chosen so that the spectra correspond to typical spectra of 
low/hard state of GBH; kTe = 100 keV, Tes = 0.8. The bulk ve- 
locity is /3bulk = 0-3. The dashed line represents the soft photon 
input, solid line shows the comptonized continuum scattered to- 
wards the observer (angle averaged), the short-dashed line shows 
the backscattered spectrum and the short-long-dashed line rep- 
resents the reflected spectrum for ionization parameter, § = 10 . 



wards the disc are still necessary in some cases in order to 
explain full range of the F-R relation found by Zdziarski et 
al. (1999), as the enhancement of flux towards the disc due 
to plasma bulk motion is required to explain R > 1 seen in 
some sources. 



3.2 Radiation spectra of outflowing corona 

In Figure H we present an example of the overall spec- 
trum calculated from the model of the outflowing corona. 
We choose the following values for the model parameters: 
k% = 100 keV, r = 0.8, kTo =0.1 keV and the bulk velocity 
/3buik = 0.3. Adopted soft photon temperature corresponds 
to a typical value for galactic black holes. Such a parameteri- 
zation is convenient if we do not consider the energy balance 
within the corona. We assume that the corona is a contin- 
uous medium, i.e. we neglect the dumpiness of the corona 
described by the parameter Hs in the model of Beloborodov 
(1999) since we also neglect the secondary reprocessing of 
the reflected component through the corona. 

The continuous line shows the radiation emitted to- 
wards an observer (averaged over the entire hemisphere, 
roughly corresponding to an inclination of 60°) while the 
short- dashed line shows the component backscattered to- 
wards the disk. Since the first scattering dominates soft X- 
ray band for galactic sources, the backscattered radiation in 
this band is enhanced, as predicted by analytical results pre- 
sented in Section 0. However, hard X-ray part is dominated 
by multiple scattered photons, anisotropy effect is smeared 
off and the backscattered component is not enhanced in this 
band. The net effect is therefore the systematic difference 
between the spectral slope of the back-scattered radiation 



and forward-scattered radiation. This effect was discussed 
for a corona without a bulk motion in a number of papers 
(e.g. Stern et al. 1995). Therefore, the continuum formed 
in the corona and emitted towards an observer is slightly 
curved, particularly in the soft X-ray band, instead of being 
a simple power law with a high energy cut-off, as frequently 
assumed in spectral analysis of the data. 

The backscattered continuum is subsequently reflected 
by the disk surface which in our calculations is assumed to 
be ionized (£ = 10*). We show this spectral component in 
the Figure |8[ 

The reflected component in the outflowing corona 
model is again partially reprocessed by the corona. The ef- 
fect depends on the corona dumpiness since only a fraction 
of radiation /is (following the notation of Beloborodov, 1999) 
would pass again through the hot plasma. In the present pa- 
per we neglect this secondary reprocessing since it is essen- 
tial only if fis is close to 1 and the optical depth is close to 
1. However, in detail modeling this effect should be rather 
taken into account. 



4 DISCUSSION 

It is generally assumed that in galactic X-ray sources and 
active galactic nuclei soft X-ray radiation originates from the 
geometrically thin and optically thick accretion disc while 
the hard flux is produced by Comptonization in optically 
thin plasma outside the disc. Part of the hard X-ray flux that 
is directed towards the observer is detected in the form of 
power law continuum and the part of the flux scattered back 
to the disc surface produces so called 'reflection hump' in the 
spectrum over 10 keV (Lightman & White 1988; Pounds et 
al. 1990; Done et al. 1992) as well as the iron Ka line near 
6.4 keV (Zycki & Czerny 1994). 

The observed values of the amplitude of the reflec- 
tion component cover broader range than initially expected 
(i? ~ 1 for AGN in Pounds et al. 1990). For some X-ray 
sources the observed reflection seems to be weak, which was 
modelled either by the disruption of the inner disc and X- 
ray irradiation of its outer parts, or by the high ionization 
state of illuminated medium, which results in steepening of 
the reflected continuum and makes it indistinguishable from 
the primary power law (Ross et al. 1999). 

Both interpretations can explain the observed correla- 
tion between R and F. In the model with a disrupted cold 
disc partially overlapping the innermost hot flow the ob- 
served correlation is due to radiative coupling of the two 
components, as the amount of overlap varies (Zdziarski et 
al. 1999). In the model with strongly ionized disc surface 
the correlation is governed by the variable optical depth of 
the ionized scattering layer (Nayakshin et al. 2000), which 
determines the effective albedo of the disc and the soft flux 
from thermalized fraction of the illuminating X-rays. Both 
models can explain the reflection amplitude values R ^ 1. 

However, the strength of reflection was for some sources 
found to be i? > 1, which may either indicate that reflecting 
medium subtends a solid angle larger than 27r, or that the 
radiation directed towards the disc is enhanced. The latter 
is possible when the scattering process is anisotropic or may 
be due to the velocity of systematic bulk motion directed 
towards the disc. 
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In this article we show that the anisotropy effect is im- 
portant mostly for the first scattering and weakens in numer- 
ical Monte Carlo simulations performed for multiple scatter- 
ing in the plasma parameterized by optical depth and elec- 
tron temperature. The maximum amplitude reached in the 
case of highly ionized reflector is only R ~ f.l. Therefore, 
the systematic mildly-relativistic motion towards the disc, or 
the 'coronal inflow', is required to explain the higher values 
of R. On the other hand, in order to produce _R > 1 in the 
disrupted disc model it would be necessary to allow for re- 
flection from outer regions of a thickened disc (e.g. Shakura 
& Sunyaev 1973) and/or from the dusty/molecular torus 
around central black holes of AGN. Since the absorption of 
irradiating X-rays below --^ 5 keV is reduced for a highly 
ionized disc surface, the reflected photons may contribute 
to the soft X-ray excesses observed in the spectra of GBH 
in the hard/low state (Cyg X-1; Di Salvo et al. 1999, and 
in preparation). We show that this contribution is higher 
when the radiation backscattered to the disc is enhanced as 
a result of the anisotropy of the first scattering. 

The hard to low state transition characteristic for many 
accreting black hole systems is in the 'disc plus sphere' model 
connected with the change of the inner radius of the cold 
disc. In the hard spectral state the cold disc is pushed out- 
ward while the hot inner plasma is responsible for hard X- 
ray radiation. In the soft state the thermal disc emission 
dominates as the cold disc extends almost to the marginally 
stable orbit. However, the physical mechanism of such be- 
haviour is currently unclear. In the outflowing corona model 
the outflow velocity is the basic control parameter. The hard 
spectral state would then correspond to rapid coronal ex- 
pansion, while the spectrum with dominating soft compo- 
nent would be produced during the vertical collapse of the 
hot gas. Again the physical mechanism of such dependence 
is unclear, and in particular it is unclear how the required 
changes of either Rin or /3buik could be driven by changing 
accretion rate. 

In conclusion, the basic predictions of all the proposed 
models are similar and only detailed computations of spec- 
tral and temporal behaviour and comparison with the high 
quality data, as expected from Chandra and XMM, may 
allow for a distinction between them. 
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APPENDIX A: THE MONTE CARLO CODE 

Our initial comptonization code was written following 
closely descriptions given by Pozdnyakov, Sobol & Sunyaev 
(1983) and Gorecki & Wilczewski (1984), and applied to 
model X-ray spectra of accretion disc with accreting advec- 
tive corona by Janiuk, Zycki & Czerny (2000). In order to 
include the effect of bulk motion of the plasma, we needed 
to make two major modifications. 

Firstly, the average scattering cross section has to be 
modified. The escape probability of a photon is given by 
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P(d)=exp( - / N,{a)dl\ , (Al) 

where A^'o = J N{v)d^v is the electron density, N{v) is the 
electron velocity distribution, d is the distance to the cloud 
boundary along the direction of photon motion, ft, and 



(a) = 



N{v){l-v-n/c)a{x)d'^v 



(A2) 



is the scattering cross section averaged over N{v). Here 



a; = 2-^7(1 - i; ■ n/c) 

TTleC^ 



(A3) 
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is the energy of an incoming photon in the electron rest 
frame, 7 is the Lorentz factor and (t{x) is the Klein-Nishina 
cross section. For an isotropic N{v) e.g. Maxwell distribu- 
tion, (a) is a function of photon energy only (for a given 

fcTe). 

With the non-zero bulk velocity we cannot use the 
method presented b y P ozdnyakov et al. (1983) to evaluate 
the integral in Eq. (A2), since N{v) is no longer isotropic. 
The presence of the specific direction - the bulk velocity 
vector /3 - introduces an additional angular dependence of 
(a) . We now have to compute the integral as in Eq. (JAj) but 
with V in the dot-product v • ft replaced by total electron 
velocity, u. Here u is the sum (in the sense of Lorentz trans- 
formation) of the thermal velocity and the bulk velocity. We 
compute this 3-D integral numerically. Introducing a coor- 
dinate system with the z-axis along the bulk velocity and 
the a;-axis along the direction of photon motion we obtain 
f2 = (sinS, O,cos&), where 6 is the angle between ft and /3. 
Applying the Lorentz transformation we obtain the electron 
velocity in the disc frame, 






(A4) 



with 7 = (1 — /S^)"""^' ^, and «x, t'y and v^ ~ components of the 
electron's thermal velocity. This enables to compute u ■ $7 
and calculate the required integral, which is now additionally 
a function of photon's direction of motion, 6. We used the 
procedure quadSd from Press et al. (1992) to evaluate the 
integral. 

The second modification concerns the scattering event. 
Since this is modeled in the electron rest frame, we intro- 
duced a pair of additional Lorentz transformations of a pho- 
ton's momentum: from the disc/corona frame to the frame 
comoving with the bulk velocity before simulating the scat- 
tering event, and the reverse transformation after the scat- 
tering event. 

This paper has been processed by the authors using the 
Blackwell Scientific Publications JM^jK style file. 
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